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In this paper we explore the incommensurate spatial modulation of spin-spin correlations as the 
intrinsic property of the doped Mott insulator, described by the t — J model. We show that such an 
incommensurability is a direct manifestation of the phase string effect introduced by doped holes in 
both one- and two-dimensional cases. The magnetic incommensurate peaks of dynamic spin suscep- 
tibility in momentum space are in agreement with the neutron-scattering measurement of cuprate 
superconductors in both position and doping dependence. In particular, this incommensurate struc- 
ture can naturally reconcile the neutron-scattering and NMR experiments of cuprates. 
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I. INTRODUCTION 

Antiferromagnetic (AF) spin correlations in cuprates persist into metallic phase as evidenced, for example, in 
neutron-scattering measurements. Compared to insulating phase, however, a qualitative change in the nature of 
AF correlations has also been observed-Efti The original magnetic peak in dynamic spin susceptibility centered at 
momentum (jr, tt) splits into four incommensurate peaks.at (tt ± ttS, tt) and (tt, tt ± tt5) in metallic or superconducting 
regime, as has been found in the LSCO compoundsEju (a recent measurements indicates the similar structure for 
YBCO compounds as well). The incommensurability are insensitive to both energy and temperature, corresponding 
to a new spatial modulation of spin-spin correlations in real space. 

The incommensurability-,,^ has been identified to be proportional to doping concentration S in approximately a 
universal fashion: 6 ~ 26 .eU| iThis observation poses a critical challenge to theories of incommensurability based on 
the Fermi-surface topology:U~El since the Fermi-surface encloses a 1 — 6 amount of electrons according to the Luttinger 
theorem, it is not a trivial issue why the doping concentration 6 of holes should dictate the incommensurability. A 
very similar issue has already been well known in the transport properties of cuprates where charge carriers are always 
characterized by hole-nature with the total number determined by 5 despite a large Fermi surface. This phenomenon 
has been interpreted due to the existence of a doped-Mott-insulatora because of the strong Coulomb interaction 
among electrons. With striking similarity, magnetic incommensurability may be another intrinsic consequence of the 
doped-Mott-insulator. 

A further challenge to any theory of magnetic-incommensurability comes from NMR spin relaxation rates, X/T\. As 
pointed out by Walstedt, Shastry, and Cheongjld if the non-Korringa behavior of 1/Ti for planar 63 Cu nuclear spins 
— a hallmark of normal-state anomalies of cuprates — originates from the AF correlations, this effect may generally 
'leak' to the planar 17 O nuclear spins if the AF correlations become incommensurate. Thus the non-Korrigna. behavior 
would have to show up in the temperature dependence of l/ 17 Ti. This was indeed predicted by theoriesou where 
the incommensurability is interpreted in terms of the FermUsurface topology. But experimentally the 'leakage' of 
non-Korringa behavior to 17 O spins has never been observed,E£l and the temperature behavior of l/ 17 Ti is drastically 
different from that of 1/ Ti, implying a near-perfect commensurate AF correlations within a one-component-theory 
framework. So how to reconcile the NMR and neutron measurements remains to be a serious puzzle in constructing 
a consistent microscopic theory. 

In this paper, we present a theoretical understanding of spatially modulated spin-spin correlations as a result of 
the doped Mott insulator described by the t — J modaL i-Such a Mott insulator at half- filling is an antiferromagnet 
whose magnetic properties have been well understood.EjE£l A hole doped into such-an antiferromagnet will induce a 
phase string 'defect' in the spin background which is nonrepairable at low-energy.E3 It represents a singular doping 
effect and profoundly modifies both spin and charge correlations in the presence of finite density of doped holes. Ej At 
low-doping where the 'rigidity' of the haltfilled AF state still persists, described by a bosonic resonating valence bond 
(RVB) pairing Jij'EEl a mean-field theoryEZl incorporating the phase string effect is able to combine the AF insulating 
phase, superconducting phase, and an anomalous normal state within a unified phase diagram naturally. In the 
present paper we demonstrate that it is in this regime a spatial modulation of spin-spin correlations indeed occurs 
as a direct manifestation of the nonrepairable phase string effect. Such an incommensurability is consistent with the 
experimental measurements in cuprate superconductors. In particular, such a mechanism for spatial modulation can 
uniquely reconcile the above-mentioned neutron and NMR experiments in cuprates. 



1 



The present work will be based on the so-called phase string formalism^ of the t — J model in which the phase 
string effect is explicitly incorporated into the Hamiltonian. In this formalism, spin raising operator reads 

s+=&: t m-i)v^. 

Here bi a denotes a bosonic spinon annihilation operator. Compared to the conventional Schwinger-boson formalism,lll 
a distinctive feature is the phase <5>' 1 which is defined by 

S^^Im In (1-2) 

It describes a phase vortex structure centered on holons: nf is the holon number operator at site I and Zi — Xi + iyi 
represents the complex coordinate of a lattice site i. If one changes z, continuously around a holon site I once, $f 
will gain an additional 2ir phase. Such a vorticity around a holon reflects the topology of the phase string effect as 
schematically illustrated in Fig. 1. In this way, spins become nonlocally dependent on the positions of holes and the 
doping effect explicitly enters spin-spin correlations through $J\ in addition to a modification of spinon spectrum via 
the Hamiltonian. 

In Sec. II, we show that a n in commensurate spatial oscillation because of will modulate the original commensu- 



rate oscillation of (— l) 1 in (1.1) at finite doping in both one-dimensional (ID) and two-dimensional (2D) cases. The 
incommensurability in ID recovers the well-known Luttinger liquid behavior, while in 2D it explains the experimental 
results in cuprates including the incommensurate peak positions in momentum space and the doping dependence of 
the incommensurability. We then further show that the effect of such a spatial modulation of spin-spin correlations 
cannot be picked up by NMR spin relaxation rates precisely due to its singular nature associated with doped holes, 
and thus reconcile the neutron and NMR experiments in cuprates. Finally, a conclusive discussion is presented in Sec. 
III. 



II. INCOMMENSURABILITY DUE TO PHASE-STRING EFFECT 



According to (1.1), the transverse spin-spin correlation function can be written down as follows 

(S+(t)Sr(0)) = (-l)'-i (^(^^(^.(4)6^(0)6^(0)) , (2.1) 

where 

Il^iJse^We-**^ ). (2.2) 

(The longitudinal correlation function will be discussed in the Appendix.) At half- filling, = and the prefactor 
[—iy-3 i n ( |2.1| ) determines a commensurate AF oscillation of the spin-spin correlation function: (— l) 4- - 7 = e l Q°'( r< ~ r " 
with the AF wavevector 

Qo = (±7r/o,±7r/o), (2.3) 

where a is the lattice constant. [Note that the average in (^J) will also produce a (small) term which oscillates 
as (— ly^i in real space aad compensates the prefactor in ( |2.1| ), leading to a ferromagnetic component of spin-spin 
correlations at half-filling 



At finite doping, the topological phase $^ defined in (1.2) will qualitatively modify the spatial modulation of the 



AF spin-spin correlations. For the purpose of illustration, in the following we discuss the ID case first as an example, 
where the phase string effect has been already showrO to play a crucial role responsible for Luttinger liquid behavior. 
Then we discuss the 2D case and make comparison with experiment. 

A. Incommensurability in ID 

In ID case, simply reduces tell 

*? = ±7r5>?. (2.4) 
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Let us consider the equal-time case first, where 



R h = e <(**-**) = gi-E^^^. (2.5) 



Here we assu me i < j without loss of generality. 



Expression (2.5) clearly shows that every holon sitting between i and j will contribute a tt phase shift to II 



Due to the spin-charge separation in ID, one may take separate averages over spinon and holon in (2.1). Define 
Nij — Ei</<j n f ■ The average of 11^- over holon degrees of freedom can be evaluated as follows 

= ^SQtx.-Xj) x g j QW osc iUation part, (2.6) 



with 



6Q( Xi - Xj ) = ±tt(N^). (2.7) 



To leading approximation, 11^ will then contribute to an additional oscillation factor at wavevector 

6Q = ±nS/a. (2. 



Besides such a singular oscillating pact, there is a slow oscillation part in (|2.6| ) which only gives rise to an additional 
decay, determined at J — > + limit ascJ ~ l/jxyl 1 ^ 2 at \xij \ = — Xj\ 3> a, and thus 



e 



iSQxi 



nfi) ~ t-^- (2-9) 



It then leads to an incommensurate momentum structure at Qq + SQ = ±(7r/a)[l — 5} = ±2kf in ID spin-spin 
correlation function. Including the time t in ( |2.9| ) only changes the decay factor \xi 3 -\" 1 / 2 to (xfj — v^t 2 )^ 1 ^ 4 at long- 
time limit (v c is the Fermi velocity of holon) without modifying the incommensurate position at 2k f. The correct 
power-law decay of the spin-spin correlation function at J/t — > + limit is given by: 

^■ w-ww-w" (2 ' 10) 

where v s is the spinon velocity. Furthermore, the similar phase-string effect also gives rise to a Luttitiger-liquid type 
Fermi-surface singularity in the momentum distribution of electrons at kf with an exponent l/8.Eil Even though 
the 2k f incommensurate structure in spin-spin correlations is still connected to the Fermi- momentum kf here, a 
Fermi-liquid type interpretation in terms of the Kohn anomaly is no longer accurate since the electron quasiparticlc 
description is not elementary anymore. A correct picture is that each holon as a topological object carries a spin 
domain- wall with it so that in calculating the large-distance spin-spin correlation function one has to count how many 



holons in between as each of them contributes a tt phase-shift as shown in (2.5). 



B. Incommensurate structure in 2D 



The spin-charge separation in 2D is slightly different from ID. While in ID spinons and holons can be treated as 
decoupled at low-energy and long- wave length limit, spinons and holons in 2D, though they still belong to independent 
degrees of freedom, will influence each other through some topological gauge fields even in_Lpng-distance limit. For 
example, spinons are subjected to a lattice gauge field A^ m ^ in 2D in the HamiltoniarO (he re (Im) refers to a 



nearest- neighbor link). As a physical quantity, the spin-spin correlation function defined in (2.1) should be gauge- 
invariant, independent of the gauge choice of A^ m . This requires Ily to be able to absorb the gauge phase arising 
from the spinon propagator Jinder a gauge transformation: Af m — > A\ L m + (8i — m ) and bi a — > 6; CT e io ' e ' which leaves 
the Hamiltonian unchanged O 

In fact, itself can be expressed in terms of a gauge field. Introduce the following relation 

/2 x (r i*0 ) 
dr ■ — j — — , (2-11) 
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where T is an arbitrary path connecting i and j without crossing the site I. In the following we always choose T as 
one of the shortest-paths between i and j without crossing any lattice site. Then fl' 1 - at t — can be expressed as 



IT 



where 



X (r- !•;) 

' |r-r,P • 



(2.12) 



(2.13) 



Similarly the lattice gauge field Af m in the Hamiltoniar 
follows 



can be also expressed by such a vector potential A h (r) as 



Ah, _ 
A lm — 



j dv-A h {r). 



(2.14) 



Then it is straightforward to verify the gauge invariance of the spin-spin correlation function defined in (2.1). 

In the mean-field theory, the gauge field Af m defined at a nearest-neighbor link is treated by smearing the holon 
number distribution locallyO It is equivalent to replacing A h (r) in (2.14) by its continuum version A' l (r) defined by 



A»(r) = \ 



(2.15) 



where the holon density p h (r') is obtained by smearing the holon-number distribution on lattice site, nf, in the 
continuum space at a scale of a. However, if one naively does the same replacement in II^- [( p. 12 )] for a larger spatial 
separation of i and j, an important effect similar to the singular phase string effect discussed in ID would be lost. 
Generally one may define 



A4n 



2 / r *- 



A h (r)- A h (r) 



and rewrite 
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h _ il f dr-A h (r)+iA0 r 
_■ — P, J r 



(2.16) 



(2.17) 



A0r will then keep track of the aforementioned singular effect in passing from A^(r) to A /l (r). 

According to (2.11), the contribution to 2 J r dr ■ A h (r) from a holon sitting at site Z is gi ven by the angle spanned 



between two straight lines connecting l-i and l—j as shown in Fig. 2a. The r.h.s. of (2.11) is a continuous function 
of r; as long as r/ stays away from the path T connecting i and j. So when I is shifted at a scale of a in the limit 
of — Tj\ ^S> a, the change of the corresponding angle is negligibly small. In this case, taking a continuum limit by 
neglecting the u nderly ing lattice can be justified. But if r/ crosses the branch-cut F, then it is easy to see that the 
angle de fined in ( 2.11 ) has a 2ir jump. Fig. 2a shows a nearest-neighbor lattice link that is across the path T and the 
value of (2.11) changes from ~ ir to —ir, namely, by ~ 2ir for the two lattice sites of the link. It is important to note 



that in II^~ such a 2ir change has no effect and thus a holon sitting at either lattice site of the nearest-neighbor link 
across T has approximately the same contribution to Hy by a phase shift ir. On the other hand, A h (v) is obtained 
by taking a continuum limit with holon distribution being 'coarse-grained' at scale of a, and in 2 L dr ■ A h (r) the 
contribution from those holons near the path T within the lattice scale of a will be simply cancelled out due to the 
opposite signs across T as shown in Fig. 2a. Thus, we find 



J2 J r dr-A h (r) ^ i2 dr- A h (r)±iTrN h (r) 



(2.18) 



where A^^(r) denotes the total number of holons at those nearest-neighboring links across the path T shown in Fig. 
2b. 

So the nonlocal phase shift in 2D is identified by 

A0 r =s ±7rN h {T). (2.19) 

As shown in Fig. 2b, the path T is always 'covered' by nearest-neighbor lattice links along the y axis if |xj — Xj\ > 
\Ui — | - Noting that each site the average holon number is S, one has 
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(N h (T)) ~ 26 ^ (2.20) 
a 

Similarly, the path T in Fig. 2b can be covered by nearest-neighbor lattice links along x axis if \xi — Xj\ < \yi — 
and the effective holon number along the path Y is estimated to be 2S\yi — yj\/a. Thus, two kinds of oscillations in 
real space are combined to give 

e^r „ 1 {g* C0S [ 27 rS/a(x t - ay)] + 0». cos[2ir6 /a(y t ~ %)]} , (2.21) 
where 9*j y ^ — 1 if \xi — Xj\ > (<)\y% — Vj\ and otherwise 0^^ = 0. 



Finally, the spin-spin correlation function (2.1) is expressed by 



(S+(t)Sr(0)) c e <Qo-(r«-*)-HA*r (bl imi(t) (>J>A^ 6 t ;(0 ) 6jT( o)^ . (2.22) 

The singular phase string effect is thus retained in the phase factor e iA ^ r and the gaug e-invariant quantity inside (...) 
on the r.h.s. can be now evaluated by the mean- field approximation. In obtaining ( 2.22 ) we have made approximations 



by neglecting the fluctuations in Atfir itself as well as the temporal dependence of life (t) . These effects in ID lead 
to an additional decay of the spin-spin correlation function in spatial and temporal space but do not change the 
incommensurate momentum position. We expect the same thing to happen in 2D where these effects merely cause 
some extra broadening of magnetic peaks. These broadening effects should become minimal atJow-temperature when 
a Bose condensation of holons occurs, which corresponds to a superconducting condensationEH (in the phase string 
theory both spinon and holon are bosonic). 

Therefore, similar to ID, the spatial oscillation factor e l A( ^ r in (2.22) will lead to in comm ensurate momentum shifts 



from Q by SQ = (±5n/a, 0) and (0, ±Sn/a) according to (2.21) [It is noted that in(2.21) will modify the plane- wave 



oscillation function but will not change the peak positions at 5Q]. The incommensurability S = 26 obtained here is 
in agreement with the experiments discussed in the beginning. Of course, whether such an incommensurate structure 
can become obser vable also depends on the broadening effect introduced by the spinon propagator [the average on 



the r.h.s. of (2.22)]. With the singular phase string effect being explicitly sorted out, such a gauge-invariant quantity 
can be directly calculated in terms of the mean-field theory developed in Ref. [l7]. For the sake of compactness, the 
detailed mat hema tical manipulation is given in the A ppendix. The dynamic spin susceptibility function x"(hji LU ) 
is shown in (A24J) (transverse component) and (A29) (z-component) and is rotationally invariant at the mean- field 



level. 

In Ref. [T^, two kinds of mean-field solution have been found for metallic phase: a uniform phase without a Bose 
condensation of spinons and an inhomogencous phase with the Bose condensation of spinons at low temperature. 



Mathematically they are controlled by the strength of fluctuations in the gauge field defined in ( A7). Note that 



the Bose condensation of spinons at half-filling corresponds to a long-range AF ordering of spins. In the metallic 
phase, the long-range AF order is generally absent. But in the spinon Bose condensation case, a tendency towards the 
AF ordering is still present as exhibited in a form of inhomogeneity (microscopic phase separation) and is argued to 
occur usually at smaller doping (underdoped regime). An experimental characterization of those two kinds of metallic 
phase has been also proposed based on the distinct energy structural!] of the local dynamic spin susceptibility function 
x'l(ui), which is illustrated in the insert of Fig. 3 for the same doping concentration 5 = 0.143 and at T = 0. The 
dashed curve corresponds to the uniform phase where at low energy there only a resonance- like peak at u> = 2E S ~ 0.4 J 
(J is the superexchange coupling constant) and a real spin gap below the peak; the solid curve corresponds to the 
inhomogeneous phase with a spinon Bose condensation where a double-peak is shown with a nonzero but suppressed 
weight at low energy resembling a pseudo-gap picture. 

We compute the momentum-dependence of the dynamic spin susceptibility function corresponding to those two 
cases with the same parameters. Fig. 3 shows x"(q, cj) as the Fourier transformation of x"(i,j;u)) given in the 
Appendix. The momentum scan is illustrated by the left top insert in Fig. 3 which is along g x -axis at a fixed 
q y = n/a, and four full circles denote the incommensurate peak positions at Qo + SQ. First of all, the dashed curve 
corresponds to aforeme ntion ed uniform- phase case at u> = 2E S . It shows no explicit incommensurate peak splitting 
because the average in (|2.22[ ) contributes to a large broadening in momentum space which smears the peak structure 
composed of four incommensurate peaks. The corresponding resonance-like peak around u> =J2Ji! s in the insert of Fig. 
3 was used to explainEj the 41 meV peak in the optimally doped Y Ba^Cu^O^ compound.Efl A single it broad AF 
peak shown in Fig. 3 is also consistent with the neutron-scattering experiments in such a compound. 

The incommensurate peaks at Qo + <5Q explicitly show up in the solid curve (inhomogeneous phase) in Fig. 3 
at small u). Since there is a spinon Bose condensation in this case, the susceptibility function has nonzero weight 
persists over to low energy until co = (the insert of Fig. 3). In this case, the spin-spin correlation length at small 
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to becomes much longer than the uniform-phase case such that the underlying incommensurate structure becomes 
visible. Experimentally, the incommensurate splitting has been indeed observed at small energy transfer in LSCO 
compoundsfina at the same positions and with the same doping-dependent incommensurability 5 shown in Fig. .-3. 
Most recently, a similar incommensurate structure has been also identified in the underdoped Y BCO compounds.El 



C. Reconciliation of the incommensurability with NMR measurements 

The NMR experiment has served as an another powerful tool in probing the magnetic properties in cuprates. A 
combination of NMR and neutron-scattering measurements should provide a more complete picture for the nature of 
AF spin-spin correlations in these materials. 

The NMR spin-lattice relaxation rate for nuclear spins at planar 17 O sites, l/ 17 Ti, can be expressed in terms of 
the real-space spin-spin correlations as followsEjEBl 



-E 



G, 



- E 

j=nn{i) 



G, 



(2.23) 



where 



Gi 



(2.24) 



According to ( p. 23 ), l/ 17 Ti only involves real-space spin-spin correlatkms-wp to nearest-neighboring sites. The spin 
relaxation rate for planar 63 Cu nuclear spin can be similarly obtamcppE3 which involves spin-spin correlations up 
to next-nearest-neighboring sites. So in the one-component theoryEjLJ of the NMR in cuprates, only short-distance 
spin-spin correlations are pjpjant. 

It has been well knownEilcj that a sharp commensur ate A F component in x"{hj) will have strong cancellation 
between the on-site and nearest-neighboring-site terms in (2.23) due to the staggered factor (— 1) I_J . So commensurate 
AF correlations, which presumably are responsible to a so-called non-Korringa temperatuce,rbehavior in l/ 63 Ti, will 
have much less significant contribution to 1/ T±, as has been discussed by many authors,Eac3 and this is consistent 
with experimental measurements where l/ 17 Ti shows a much weaker and quite different behavior usually attributed 
to the non-AF correlations in \" ■ 

However, the neutron-scattering measurements indicate AF correlations in cuprates should be incommensurate 
in nature as discussed in the previous section. Once the sharjs— commensurate AF component is replaced by an 
incommensurate AF structure experimentally, it was pointed outt2l that the non-Korringa behavior of 1/Tj at 63 Cu 
would usually 'leak' to the 1/ 17 sites, causing a temperature behavior in l/ 17 Ti which is not consistent with 
experimental results. This is due to an imperfect cancellation of the AF component of x" m (2.23) after the change 
from a sharp coinmensurate to an incommensurate (or a more broadened) magnetic structure around Qo. Indeed, 
various theoriesEm based on the Fermi surface topology in which the incommensurability arises from the Kohn 
anomaly have predicted a similar non-Korringa behavior for l/ 17 Xi, in contrast to the experimental results. In this 
sense, incommensurate AF correlations seem inconsistent with the NMR data. In order to resolve this issue, some 
authorstJ even have proposed to modify the form of hyperfine coupling as an alternative way out. 

But such a paradox can be easily reconciled in the phase string description without changing the hyperfine inter- 
actions which-are determined based on the local chemical structure of cuprates and are consistent with the Knight 
shi ft da ta.E^rE3 In the present theory, the incommensurate structure or spatial modulation, as represented by e lA ^ r 
in ( 2.22Q , has been directly attributed to the phase vortex structure whose singular centers are bound to holons. The 
resulting incommensurate oscillation in spin-spin correlation function is only meaningful at a scale much larger than 
the size of the singular center. As for a nearest-neighbor lattice link (ij), Eft will simply reduce to 



(2.25) 



where is the lattice gauge field given in ( 2.14 ). In the mean-field approximation is replaced by J? dr ■ A' l (r) 
as discussed before, without causing any additional singular phase shift, namely, A</>r = 0. It is because th ere ca n be 
no holon sitting between two nearest-neighboring spins. Trivially, one also has Eft, = 1. Thus, l/ 17 Ti in (g~23|) will 
simply not be able to 'see' any incommensurate phase shift, even though such a incommensurate spatial modulation 
does exist in large-distance correlations as given by (2.21). The latter can be picked up by the neutron-scattering 
measurement as discussed in previous section. 
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Fig. 4 shows an example of a non-Korringa behavior of l/ 63 Ti vs. a Korringa-like behavior of l/ 17 Ti based on 
the dynamic spin susceptibility function determined in the Appendix where the uniform (optimal-doping) mean-field 
solutioriiH is used at 6 = 0.143. The distinct temperature-dependence and the relative size of l/ 63 Ti and 1/ 17 7\ 
here are in accord with the typical experimental data of cuprates. l/ 63 TiT exhibits a non-Korringa 1/T behavior 
at high temperature while l/ 17 TiT remains roughly a constant (Korringa behavior). Note that the decrease of the 
relaxation rates at low-T is-due to the spin gap ~ E s which also decides the superconducting transition temperature 
T c around the same scaleJHI One thing we need to emphasize here is that even though in the uniform-phase case 
the momentum structure in Fig. 3 looks like a commensurate peak, it actually is composed of four incommensurate 
peaks with stronger broadening, and thus the 'leakage' effect would be expected even stronger in this case. In fact 
if we assume the incommensurate spatial modulation ( [2.21 ) to persist into short-distance ~ a, the non-Korringa 
temperature dependence will then immediately show up in l/ 17 Ti precise ly a s pointed out in Ref. |TT^. Therefore, the 
origin of the incommensurability from the nonlocal vortex structure of (1.2), i.e., the phase string effect, is indeed 
crucial to resolve this issue. 



III. DISCUSSIONS 



In the present paper, we have related the magnetic incommensurability to the intrinsic properties of the doped Mott 
insulator. In the Mott insulator, the spins are the only degrees of freedom not frozen by the strong on-site Coulomb 
interaction. Holes doped into such a Mott insulator always introduces a phase frustration (phase string effect) on 
the spin background by violating the Marshall sign rule.E3o Such a phase string effect then contributes a unique 
spatial modulation to the spin-spin correlations as demonstrated in this paper, which is more or less independent 
of energy and temperature but linearly depends on the doping concentration. In momentum space it leads to an 
explicit incommensurate splitting of the AF peak in the dynamic spin susceptibility when the spin correlation length 
is sufficiently long, in agreement with the neutron-scattering measurements of cuprates in both the peak positions and 
doping dependence. Furthermore, due to the phase string origin of this incommensurability, it has been shown that 
the NMR spin relaxation rates should not be able to pick up the incommensurability AF correlations at length scales 
comparable to the lattice constant, and thus conflicting implications of neutron and NMR data can be reconciled 
within the present framework. Conversely, one may also conclude that the magnetic incommensurability observed 
in cuprate superconductors provides a direct experimental evidence for the phase string effect in the doped Mott 
insulator. 

We conclude this paper by making several remarks on some related issues. First of all, we point out that the present 
magnetic incommensurability is a unique property of metallic phase, which should disappear in the localized regime. 
Recall that the phase string effect leads to a vortex phase e 1 ** emerging in spin operator (1.1) which is the reason 
for the incommensurate spatial modulation. It has been previously discussecO that the same phase e I$i also causes 
the disappearance of the AF long range order. It was argued that only in the hole-localized regime (insulating phase) 
where the phase string effect is no longer effective, the phase e I$i can be compensated by the vortex phase carried by 
spinons which become vortices in the insulating phase, and the AF long range order may be recovered. By the same 
token, the incommensurability induced by e I$i should be also gone in this localized regime at small doping. This is 
in accord with the experiments where only a commensurate AF peak has been observed in the insulating phase.D 

Secondly, as briefly mentioned in Sec. IIA, even though 2k f incommensurate oscillation of ID spin-spin correlation 
function may be viewed as the Kohn anomaly of the Fermi surface effect, it is not very meaningful because single 
electron excitations are no longer elementary and the spin-spin correlation function cannot be reasonably approximated 
by a 'bubble' diagram composed of the single-electron Green's function without including high order corrections. 
Similarly the 2D magnetic incommensurability in the phase string formalism does not necessarily exclude the relevance 
of the Fermi surface topology. Instead it implies that the latter explanation is no longer accurate as holon and spinon 
become new elementary excitations. Especially, the interpretation of the doping-dependent incommensurability and 
the reconciliation of neutron and NMR measurements can only be naturally realized in this phase string description 
of the doped-Mott-insulator. 

Finally, we emphasize that the incommensurability discussed in the present paper is due to the topology associated 
with each individual holon. It has nothing directly to do with the phase separation oiuiharge inhomogeneity at a larger 
length scale involving many holes as discussed in literature (like the striae phasecj), including the inhomogeneous 
(undcrdoping) metallic regime studied within the phase string formalism JI3 Nevertheless, there is still a fundamental 
connection: all of these may be attributed to the competition between the kinetic energy of holes and superexchangc 
energy of spins and the resulting tendency towards phase separation. The vorticity surrounding each holon may be 
regarded as a microscopic version of phase separation against a uniform spiral phase.l 
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APPENDIX: DYNAMIC SPIN SUSCEPTIBILITY FUNCTION 



1. Mean-Field Theory 

The t — J model in the phase string formalism can be written in two terms in the mean-field approximation:0 
H a + Hh- H s governs the spinon degrees of freedom and Hh determines the holon degrees of freedom. In the present 
paper, we mainly focus on the spin magnetic properties decided by H s : 

H s = ^2 E m ~fl ncr j m(T + constant, (Al) 

mo 

where 7^ is the creation operator of bosonic spinon elementary excitations. The spinon excitation spectrum E m in 
this mean-field theory is given by 



E m = V^-C (A2) 

in which, A m is defined by 

A m = A-a|£ m |, (A3) 

where A is the Lagrangian multiplier to be determined by the condition of total spin number equal to N(l — S) (N 
is the total number of lattice sites.) The coefficient a is a ratio of two energy scales: a = Jh/ Js, in which Jh = St' 
describes an effective hopping effect on spinon degrees of freedom with t' ~ J and J s — (J/2)A S with A s being the 
bosonic RVB order parameter characterizing the present mean-field state. The Lagrangian multiplier A and the RVB 
order parameter A s are determined by the following mean-field equations, respectively, 



ly^coth^i, (A4) 



2-5 

N^E m 2 
and 

A s = -i-V JiLcoth^. (A5) 

an ^ j s E m 2 v ' 

m 

The spectrum £ m will be modified by doping effect through a gauge field A« as the eigenvalue of the following 
equatiorO 



Uw m *(i) = -Js e toi4 «iiwO')- ( A6 ) 

Here A^- (with j = nn{i)) is originally defined by the following gauge invariant condition 

E4=^ (A7) 

c 

where C is an arbitrary counterclockwise closed path and Nq is the number of holes enclosed by it. In the present 
mean-field theory, the dynamic effect of is replaced by a randomness as it comes from an independent degrees of 
freedom (holons). Nq may be rewritten as N^ + SNq with fluctuating 8Nq. Correspondingly describes a uniform 
fictitious magnetic field with a random flux fluctuation on the top of it. Two kinds of mean-field solutions due to 
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different strengths of the fluctutations are found,0 which are distinguished by whether there is a Bose condensation 
of spinons or not at low temperature. The former case is called inhomogeneous metallic phase and the latter is called 
uniform metallic phase in Ref. 

Finally, the Bogoliubov transformation relates 7-operator to the spinon operator b{ a is given as follows: 

bia = ^2 ( U "l7mtr - wi-J ^° Xm W ma (*) , (A8) 

m 

where u m — 1 / \/2(X m / E m + l) 1 / 2 and u m = l/\/2(X m /E m — l) 1 / 2 sgn(^ m ). The phase factor e lcrXm is determined up 
to a change 

e **» _> - S gn(Z m ) x (A9) 
each time when a holon changes sublattices.0 



2. Dynamic spin susceptibility 

Let us first consider the transverse component defined in Matsubara representation as follows 

AS 

X+-(i,j;iwn)= / iT ^{^(^(0)}, (A10) 
Jo 

Here w„ = and 

(T T 5+(r)5-(0)) = (-!)<-> (r^^Tj^frjn* (r)6} ; (0)6 JT (0)) . (All) 



According to discussions in Sec. IIB, the r.h.s. of ( |All|) can be expressed by 

(-ir j (r^^^^^fet^O^^e^/r^W+^y (A12) 
In the mean-field approximation one has 

V T ^ T (r)6 li (r)6 j t i ( )& jT (0)) s - (T T 6j T (r)6j ; (0)) a (T T 6 iT (0)6 ii (7 

+ (r T 6j T (r)6 iT (0)) a (T T 6 a (r)6 3 t i (0)) s , (A13) 

where (...) s denotes the average over the spinon degrees of freedom. According to the mean-field theory outlined in 
the above section, one has 

bl(T)=e H *T b l e -»,r 

= E ( u ™7Le EmT ~ v m lm-„e- E ^) w^(i)e-^ . (A14) 

m 

Then by noting w^ n _ <7 — w ma one finds 
[T T bl{r)b\_M) 

= -^u m v m (w* ma (i)w m(T (j))e- i " A x™ (e £ '" T (7L7ma) s + e- £ '" T ( 7m -.7L- CT ) s ) , (A15) 

TCI 

at r > where e~ lljAXm is introduced to incorporate the phase shift effect ([A|). e _MTAXm may be determined based 
on how many exchanges between the spinon and background holons. The spinon will acquire a phase shift ±7r for 
£ TO > due to —sgn(t; m ) — —1 each time it exchanges positions with a holon. Thus, if denotes the number of 
holons encountered and exchanged with the involved spinon on the shortest-path T connecting i and j, the total phase 
shift acquired by the spinon at £ TO > is ztirN^ ~ A(j>r- Then we may write down 
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-iaAx„ 



1, if^ m <0 
if£ m >0 



(A16) 



It is noted that for each m with £ m < 0, one always can find a state m with £, m = — £ m > with a wavcfunction 



UW(0 = (-l) l W m(T (i) 



(A17) 



according to (A6). Then in terms of ( A16 ) and ( A17 ) we finally have (at r > 0) 

[T T bl(r)b]_A0)) s =((-ir^-l] 
x ^ u mV m {w* ma (i)w ma {j)) [n(E m )e EmT + (1 + n(E m ))e~ EmT ] , 

where the summation only runs over those states with £ m < 0. Similarly, one can also determine 

/ T T bl(r)b J A0)) s = ((-1)'-^ + l) ^ (< CT (i)«w(j)) 

m 

x [z4n(£ m )e B ™ T + + n{E m ))e- E ™ T ] , 



(A18) 



as well as other averages in (A13). 



Then after integrate out t in (A10) and identify (— l) N v with e lA ^ T in operator form, one obtains 

x+- (i, j ; iu n ) ~ xV- (h j ; + (-iy- j e iA ** ■ x ( +l (i, j ; iwn), 



(A19) 



(A20) 



where 



><(C lm 2 (l+™(£rn)+n(iw)) 1 



mm' 

l 



± , 2 n(^ m /) -n(£, 

V-r mm' ) 



iu> n + E m — E., 

1 



2 V iu> n + E m + E m ' iw n — E m — E„ 



(A21) 



with 



K + - j) = 2 

mm' V ' J ) 



E 



i \ dr-A h * / -\ / ■ \ 

e Jr W m] {l)w m ]{j) 

ia I dr-A. h * / .\ / .\ 



ia I d-A h * / ■ \ / ■ \ 



(A22) 



In obtaining the last line, the symmetry w* lcr — w m ~a has been used. The coherent factors, p mm i and l mm , are 
defined by 



3 ± 

mm 
mm 



^>,, ™ / — ^m^m/ i ^m^m' 



(A23) 



Finally, the dynamic spin susceptibility function x+— (*>J> w ) can be obtained as the imaginary part of xh after an 

analytic continuation io; n — > w + i0 + is made: 



where 



(A24) 



- E m ~ E nl ') + 2 (n(Em) - n{E m ')) {Pmm'f 5 ^ + E ™ ~ E m ,)] 

The longitudinal spin susceptibility function can be similarly obtained: 



(A25) 
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X„(i,j;iw„) = / dre^ (T T S?(r)S?(0)) 
Jo 

^xV^r,^n) + {-ir°e^ -x&frxiun), (A26) 

where Sf = a b\ a bi a IS use d and 

>n(E m >) - n{E rn ) 



(P ± 



2 \ iw n + E m + E m > iu> n — E m — E n 
with 



(A27) 



(A28) 



where the phase factors e ±I<T ^r dr A are introduced to explicitly show the gauge invariance of K^ m ,. Generally, one 

expects the gauge-invariant quantity e ia ^ drA w* na (i)w ma (j) to be independent of a or to be real, which in turn 
means K^ m , = K^^,. Consequently, the rotational invariance of spin-spin correlations is retained. The corresponding 
dynamic spin susceptibility function is given by 

X^(i I i; W ) = $y(z,i;^) + (-l) i ^ e iA ^$W(i,j;a;) ) (A29) 

where 

= \Y^ K ™m'ihj) [(l + n(E m ) + n(E m ,)) {l± m ,fsgn{u) 

min' 

■5{\u\ -E m - E m ,) + 2 (n(E m ) - n(E m ,)) {p± m ,) 2 5{uj + E m - E m ,)] . (A30) 
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Fig. 1 The topology of a holon due to the phase string: any closed loop around it will cut through the phase string 
once or an odd number of times. 



Fig. 2a The illustration of the angle defined by (2.11) at site I. Such an angle jumps approximately 27r at a 



nearest-neighbor lattice link (bold one) across the path T connecting sites i and j in the limit Ir^ — r,-| ^> a. 

Fig. 2b The path T connecting i and j is covered by the nearest-neighbor links along the y-axis (bold links) in the 
case that \xi — Xj\ > \yi — yj |. 

Fig. 3 Dynamic spin susceptibility function x"(q, w) versus q x (with a fixed q y at 7r/a) at 8 = 0.143. Two cases: 
the dashed curve corresponds to the uniform phase at ui = 2E S ~ 0.4J where xK^) ^ s peaked (the insert); the solid 
curve corresponds to the inhomogeneous phase at u> ~ 0. The w-dependence of Xl( w ) m the insert is from Ref. [Tt]. 
The positions of the incommensurability is illustrated by four full circles at the left top. Note that the maxima of \" 
in two cases are scaled to the same value for comparison. 

Fig. 4 The contrast of the non-Korringa and Korringa behavior of the NMR spin relaxation rates for planar copper 
and oxygen nuclear spins for the uniform phase case at 6 — 0.143. 
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